Pade- Improvement of QCD Running Coupling Constants, Running 
Masses, Higgs Decay Rates, and Scalar Channel Sum Rules 



V. Elias*and T. G. Steele 
Department of Physics and Engineering Physics 
^ University of Saskatchewan 

Q>^ ! Saskatoon, Saskatchewan S7N 5C6 

0\ \ Canada 

and 

F. Chishtie, R. Migneron and K. Sprague 
Department of Apphed Mathematics 
I The University of Western Ontario 

London, Ontario N6A 5B7 
. Canada 
> . 

^ ; February 1, 2008 

cn 
o 

00 ■ PACS: ll.lO.Hi, 11.55.Hx, 12.38.Aw, 14.80.Bn 

' Abstract 

c^h: . 

We discuss Pade-improvement of known four-loop order results based upon an asymptotic three- 
parameter error formula for Pade-approximants. We derive an explicit formula estimating the next- 
fH [ order coefficient R4 from the previous coefficients in a series 1 -I- Rix + R2X^ -f R^x^ . We show that 

such an estimate is within 0.18% of the known five- loop order term in the 0(1) /3-function, and within 
10% of the known five-loop term in the 0(1) anomalous mass-dimension function 7™ (5). We apply 
■ the same formula to generate a [2|2] Pade-summation of the QCD /^-function and anomalous mass 

^ , dimension in order to demonstrate both the relative insensitivity of the evolution of ctsifJ-) and the 

running quark masses to higher order corrections, as well as a somewhat increased compatibility of 
the present empirical range for as{mr) with the range anticipated via evolution from the present 
empirical range for as{Mz). For 3 < < 6 we demonstrate that positive zeros of any [2|2] Pade- 
summation estimate of the all-orders function which incorporates known two-, three-, and four- loop 
contributions necessarily correspond to ultraviolet fixed points, regardless of the unknown five-loop 
term. Pade-improvement of higher-order perturbative expressions is presented for the decay rates of 
the Higgs into two gluons and into a bb pair, and is used to show the relative insensitivity of these 
rates to higher order effects. However, Pade-improvement of the purely-perturbative component of 
scalar/pseudoscalar current correlation functions is indicative of large theoretical uncertainties in QCD 
sum rules for these channels, particularly if the continuum-threshold parameter sq is near 1 GeV'^. 



Oh! 



*Pemianent Address; Department of Applied Mathematics, The University of Western Ontario, London, Ontario N6A 
5B7, CANADA 



1 



1. Introduction 



A recent body of work [1-3] has demonstrated how higher order terms in a number of field-theoretical 
perturbative series can be estimated by Pade-approximant techniques. 

Of particular interest are applications to QCD quantities, particularly /3- and 7-functions now known 
to four- loop order in Ug [4-6]. Padc-approximant methods have already addressed the ?^/ - (flavour- 
number-) dependence of higher order terms in the QCD (5- and 7-functions [2,3]. Near cancellations 
between coefficients of successive powers of nj, however, can lead to large uncertainty in the estimated 
overall size of such higher-loop contributions - small errors in fitted coefficients of have been seen to 
lead to much larger errors in the aggregate (now known) four-loop contribution to the /3- function [2]. 

In the present paper, our focus will be on using Pade-approximant methods to estimate the magnitude 
of higher-order corrections to quantities already calculated to three- and four- loop order in QCD. We 
assess the theoretical uncertainty of such calculations by seeing how closely they coincide with their 
own Pade- improvements, as well as whether successive orders of perturbation theory exhibit convergence 
toward Pade-summation estimates of the full perturbative series. 

The present paper is phenomenologically oriented, specifically aimed at developing Pade-improved 
estimates of what are hoped to be computationally and/or experimentally accessible quantities. The 
particular items of interest considered are the 0{N) /3-function and anomalous mass dimension (Section 
2), the running QCD coupling constant (Section 3), the running quark masses (Section 4), the Higgs- 
boson decay rates into two gluons and into a bb pair (Section 5), and the purely perturbative content of 
QCD-sum rules based upon scalar and pseudoscalar current correlation functions (Section 6). In Section 
3, we also discuss some general implications of Pade-summation as an approximation to all orders of 
perturbation theory. In particular, we analyze the fixed point structure of the most general [2 | 2] 
Pade-summation estimates of the full content of QCD /3-functions for Uf = {3,4,5,6} whose Maclaurin 
expansions coincide with presently known perturbative contributions. 

In the section that immediately follows, we discuss how Pade-approximant methods can be used 
to estimate higher-order corrections to perturbative series in which the leading four terms are known. 
Although this methodology appears in other work [3], the presentation leading to eq. (2.12) [which has 
not appeared elsewhere] will, it is hoped, be of some value to those unfamiliar with Pade-approximant 
methods. We also demonstrate the reasonable agreement of predictions based on (2.12) with now-known 
five-loop terms in 0{N) (3- and 7-functions. 



2. The APAP Algorithm: Pade-Improvement for Pedestrians 

We consider the general problem of developing a Pade-improvement of the series 

S = 1 + Rix + R2x'^ + R3X^ + (2.1) 

where i22, -R3} are known and {i?4, i?5, ...} are not known, through use of the asymptotic error 
formula for estimating Rn^m+i via the Pade approximant 



_ 1 -I- aix + a2X^ + ...qnx^ 

■^[^1""! = l + biX + b2X^ + ... + bMX^ 

= 1 + Rix + R2X^ + R3X^ + ... + Rn+m+ix'^+'^'+' + ... (2.2) 
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Let R^"f^_^_i be the prediction one would obtain from the [A?^|M] Pade Approximant, and let i?jv+M+i 
be the true value of the coefficient. The structure of the asymptotic error formula is given by [2] 

R^%+,-Rn+m+i _ M\A^ 

Rn+m+1 [N + M + aM + h]^ ^ ' ' 

with numbers {A, a, b} (independent of N, M) to be determined. 

This error formula simplifies considerably if M is always chosen to be 1 [3]: the right hand side of 
(C) becomes —A/[N + 1 + (a + b)], with only the two numbers {A, a + b} to be determined. These can 
be determined explicitly for the series S in (2.1) given knowledge of the three coefficients {Ri, R2, -Raj- 
Given knowledge of Ri only, the [0|1] Pade approximant 

^fom = , \ =l-bix + bjx\.. = 1 + Rix + R^^^x^ (2.4) 
^ ' ' 1 + bix 

predicts i?!'"'^'^ = -^i- Consequently, we see from the asymptotic error formula (2.3) that 

T>Pade n 7j2 r? A 

R, -R^ ^ R^-Ra ^ -A ^ 

R2 R2 l + {a + b) ' ^ ^ 

Given knowledge of only Ri and R2, the Pade approximant 

5'[i|i] = ^^^^ = 1 + (ai - 6i)a; + 61(61 - ai)a;^ + 6i(ai - 6i)x^ + ... 

= 1 + Rix + R2x'^ + R^'''^'^x^... (2.6) 
predicts R^"'^'^ = R^/Ri- The asymptotic error formula (2.3) for this case implies that 

R2/R1 ~ Rs _ —A _ 



R3 2 + {a + b) 



(2.7) 



Given knowledge of {Ri, R2, R3} in the scries (2.1), the relative errors 62 and 83 are specified com- 
pletely by the left-hand sides of (2.5) and (2.7). These two equations may be regarded as two equations 
in the two unknowns A and (a -|- b) characterizing the asymptotic error formula (2.3) when M = 1. The 
solution to these two equations is 

A =[1/62 -1/63]-', (2.8) 

(« + 6) = |^. (2.9) 
03 - 62 

This information is sufficient to generate an asymptotic Pade approximant (APAP) estimate [2,3] of the 
unknown coefficient R4 in the series (2.1). Consider the [2|1] Pade approximant 



'[2|1] 



1 -|- aix + a2X 
~ l + bix 

= 1 -I- (ai - bi)x + [02 - 6i(ai - bi)]x'^ 
+ [-6i[a2-6i(ai-6i)]]x3 
+ [bl[a2 - b,{a, - b,)]]x^ + ... 

= 1 + Rix + R2x'^ + Rsx^ + Rf'^^x'^ (2.10) 
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The three known values of {Ri, R2, R3} completely determine the three parameters {01,02,61} charac- 
terizing >S'[2|i]- Wc sec from (2.10) that R^""*^ = R^/R2. However, the asymptotic error formula (2.3) 
suggests that a more accurate estimate of the true value R4 differs from R^"-'^^ by a predictable relative 
error: 

R, = ^'=3 + {a + by ^^-^^^ 

in which case we find from (2.11), (2.8) and (2.9) that 

R3/R2 



R4 = 



l + <54 

R2{S2 - 26i - 5162) 
R^{R2 + R1R2R3 ~ 2-R1-R3) 
-R2(2-f22 — R1R3 — -^^f-^i) 



(2.12) 



As an example, we test the applicability of (2.12) by comparing its prediction to the known 0{g^] 
coefficient of the 0(1) /3-f unction [7]: 



/3(i)(y) = 1.5^2 h _ (17/9)5 + 10.8499C/2 

- 90.53535^ + 949.523/ + 0(/)] . (2.13) 

Identifying Ri = -17/9, R2 = 10.8499 and R3 = -90.5353, we find from (2.12) that R4 = 947.8 in 
startlingly close agreement to the next {g'^) term within (2.13). Although APAP improvement has also 
been applied elsewhere [2] to 0{g^) terms in 0{N) /3-functions, the result obtained here relies on direct 
and explicit use of the full asymptotic error formula (2.3). A comparison of (2.12) predictions and exact 
values of P^, the 0{N) /^-function coefficient of g^, is presented in Table I. We emphasize that these 
predictions are not obtained by a fitting of the A?^-dependence or any knowledge of the A?^^-dependence 
of as is the case in Table 3 of ref. [2]; for comparative purposes, the predictions of ref. [2] are also 
listed in Table I. 

We can use (2.12) to predict the known R4 coefficient within the 0(1) theory's anomalous mass 
dimension [7], as well: 



7m(5) = (5/2) [1 - 0.83335 + 3.5OO5' 
- 19.965^ + 150.8/ + 0{g^ 



(2.14) 



Eq. (2.12) predicts R4 = 135.1, a result only 10% off the 150.8 value given in (2.14). Table I shows that 
predicted R4 coeflficients for ^mig) within 0(2), 0(3), and 0(4) also remain within 20% of their true 

values, as given in [7]. These results provide a reasonable basis for applying (2.12) [and its concomitant 
asymptotic error formula (2.3)] to the (3- and 7-functions of QCD, as we will do in Sections 3 and 4. 

A final improvement of the series (2.1) is possible by expressing this series as a [2|2] diagonal ap- 
proximant — this is a more accurate representation of the infinite series S than one would obtain by 
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arbitrarily truncating the series after the i?4a;^ term. Given known values of {Ri, R2, R3} and using the 
APAP estimate (2.12) for R4, the approximant S'[2|2] of the infinite series S is fully determined: 



1 + aix + a2X^ 



(2.15) 



1 + bix + b2x'^ ' 



bi = 



R1R4 ~ R2R3 

i?2 — R1R3 



(2.16) 



i?3 — R2R4: 
R2 ~ R1R3 



(2.17) 



(2.18) 



02 = -R2 + biRi + 62 ■ 



(2.19) 



Eq. (2.15), as determined from (2.12) and (2.16-19), constitutes the procedure we denote as "Pade- 
summation" of the series S in (2.1). 

3. Pade-Improvement of the QCD Coupling 

The QCD minimal subtraction (MS or MS) renormalization-group functions (3{x) and 7(2;) are now 
known to 4-loop order [4,5,6]. Prior work involving Pade- improvement methods has attempted to predict 
the flavour-dependence of these functions to 4- and 5- loop order [2,3]. Ps and /34, the 4- and 5- loop 
order corrections to the /3-function, are respectively third- and fourth-degree polynomials in n/, and Pade 
methods have already shown some success in predicting the polynomial coefficients now known for (3^ 
[2]. However, an accurate determination of the polynomial coefficients within /33 is not reflected in the 
accuracy with which is itself determined. Thus the overall Pade- driven estimate of presented 
in [2] for n/ = 3 [/?3 = (7.6 it 0.1) • 10^/256 = 30 it 1, using normalization conventions appropriate 
for (3.1) below] is substantially below the true value (/?3 = 47.23), even after allowances are made for 
claimed uncertainties arising from quadratic Casimir contributions [3]. It should also be noted that this 
estimate, once disentangled from a fitting procedure aimed at ascertaining the explicit n/ dependence 
of Ps, follows from a simplified version of the asymptotic error formula (2.3), in which the denominator 
{N + M + aM + h)^ is taken to be [2]. Indeed, it is difficult to tell at this stage whether the 

discrepancy between Pade-estimates of P^ and the true value arises primarily from quadratic-Casimir 
contributions not occurring in lower orders, or alternatively, arises from the error involved in simplifying 
the asymptotic error formula in order to make an estimate of P^ possible. Without such a simplification, 
the error formula (2.3) has (in principle) three arbitrary constants (A,a,b) instead of one (A). 

Consequently, in this section we will predict /34 directly using the APAP-algorithm (2.12) following 
from the full asymptotic error formula (2.3). We will not prejudice these predictions of Pi by attempting a 
fit of the polynomial dependence on n/, nor will we attempt to disentangle quadratic-and-higher Casimir 
contributions from P^ and /?4. Rather, we will make distinct predictions of Pi via (2.12) for = 3, 4, 5, 6. 
The validity of such an approach, particularly the possibility that the full asymptotic error formula is 
inclusive of higher-order Casimir contributions to P{x), would be best established by comparison to an 
exact calculation of /?4, when available. 

We define the /^-function as in [4] : 
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d °° 

fjf — 2^ ~ ^ Pi^ 
^ 1=0 

oo 

= -l3QX^^Rix\ (3.1a) 

i=0 

with x = as {ij) I TT and 

Ri = Pi/Po, (3.1b) 
Known values of f3o — (3^ are as follows [4,5]: 

/3o = (11 - 2n//3)/4, (3.2a) 

^1 = (102 - 38n//3)/16, (3.2b) 

(32 = (2857/2 - 5033n //18 + 325n^ /54) /64, (3.2c) 

/33 = 114.23033 -27.133944n/ + 1.5823791n^ 

+ 5.85669582 • IQ-^n^. (3.2d) 

We use (3.2), (3.1b), and (2.12) to predict the following values for (d^: 

nf = 3: i?4 = -849.74, /?4 = -1911.9; (3.3) 

n/ = 4 : i?4 = 40.203, = 83.7563; (3.4) 

nf = 5: i?4 = 70.203, /?4 = 134.56; (3.5) 

nf = 6: i?4 = -239.22, /?4 = -418.64. (3.6) 

The large negative values for nj = 3 and nj = 6 reflect the near cancellation of the factor (2i?2 ~ 
i?fi?3 — -R1-R2) in the denominator of (2.12). Since a change of sign can easily occur if this cancellation 
is over- or under-estimated, the safest interpretation of (3.3) and (3.6) is to predict a relatively large 
magnitude for /34, with the sign uncertain. 

Corresponding results from Table III of reference [3] with quadratic Casimir contributions and a 
1/1024 normalization factor appropriate to (3.1a) included are (34 = {278(n/ = 3),202(n/ = 4), 165(nj = 
5), 166(n/ = 6)}. As stated earlier, these latter results are based upon a fit to the polynomial coefficients 
of in (34, with near cancellation of very large opposite-sign coefficients for the k = and k = 1 terms. 
These latter results appear most consistent with those obtained above when nj = 5. 

To get a feeling of the magnitude of these Pade estimates of 5-loop effects, we generate the [2|2] 
approximant (2.15) from the known values for i?i,i?2, and R3 and our estimate of R4, and incorporate 
this approximant directly into the /3-function [/j^dx/dfi'^ = P{x)]: 



9x^ 

n/ = 3 : P{x) = — — 



1 + 94.383X - 75.605a;2 
1 + 92.606X - 244.71x2 



(3.7a) 
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71/ = 4: P{x) 



Hf 



~[2 



12 



1 - 5.8963X - 4.Q110x^ 
1 - 7.4363a; + 4.3932a;2 

1 - 5.9761X - 6.9861x2 
1 - 7.2369X - 0.66390x2 



(3.7b) 



(3.7c) 



We can use these Pade approximations to the full /?- function to evolve as{p,) = ttx(ij,) down to = 1 
GeV from an initial condition as{Mz) = 0.118 [8] through use of the following (4- and 5-) flavour threshold 
matching conditions with m{fj,th) = fJ-th [9]- 



th) 



th) 



1 + 0.1528 (x("/)(M)^ 



+ 



{0.9721 - 0.0847(n/ - 1)} (x("/)(Mtft))' 



(3.8) 



In Figure 1 we display the evolution of q;s(//) to low energies through use of two- loop, three- loop, 
four-loop, and the Pade-improved four-loop /3-functions of (3.7). All of these curves arc generated from 
the initial condition as{Mz) = 0.118, with /xj/j = 4.3 GeV {i.e. mh{mh) = 4.3 GeV) identified as the 
Uf = 5 flavour threshold, and with nth = 1-3 GeV identified as the n/ = 4 flavour threshold. Eq. (3.8) is 
utilized in full in both the four-loop and Pade improved calculations; it is utilized to 0{x'^) to generate 
flavour-threshold initial conditions in the thrcc-loop calculation, and the matching condition is trivial for 
the two-loop calculation. It is evident from the figure that curves from successive orders of the /3-function 
appear to converge from below to that generated via (3.7), the Pade-summation approximating all orders. 
The gaps between curves of successive order clearly narrow as the order increases. Figure 1 shows that the 
Pade-summation leads to a curve for as{p-) that exceeds the unimproved four-loop curve by less than 1%. 
Such a difference is inconsequential compared to the estimated uncertainties in as{Mz) and /ifh for four- 
and five-flavour thresholds. Both the four-loop and the Pade-improved curve are indicative of benchmark 
values as{lGeV) = 0.48 and asim-r) = 0.32. However, if the flavour thresholds and as(M^) are assigned 
their accepted [8] lower-bound values [^^(M^) = 0.115, /ijj^^'"^'' = 4.1 GeV, 41^^^^^ = 1.0 GeV], the 
values we obtain at = 1 GeV and /j, = rrir from either four-loop or Pade-summation /3-functions are 
as{lGeV) = 0.41 and 0^(772,-) = 0.29. Corresponding upper-bound values [as{Mz) = 0.121, /i^"^ = 4.5 

GeV, fj,^^-^ = 1.6 GeV] lead to as{lGeV) = 0.57 and as{mr) = 0.35. In view of these (much-) larger- 
than-1% uncertainties in low-energy values, the best possible test at present of Pade-improvement would 
be nonempirical, i.e. a comparison to an explicit 5-loop calculation of the /3-function. 

We note, however, that higher-order effects do appear to increase the overlap between the present 
(somewhat large) empirical range for as{mr) (0.370 ± 0.033 [8]) and the range predicted via evolution 
down from the present empirical range for as(Mz) (0.118 ib 0.003 [8]). Taking into account the present 
uncertainty in the five-flavour threshold {iJ,th = 4.3 it 0.3 GeV [8]) and incorporating the matching 
condition (3.8) for as{fi) below and above /i^/^, we find the following predicted ranges for 05(771,-) for 
two-loop, three-loop, four-loop and Pade-summation /3-functions: 
2-Loop: 



0.2910 < asirur) < 0.3391, 



(3.9a) 



[0:5(777^)]^ = 0.3137, 



(3.9b) 
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3- Loop: 

0.2944 < asirur) < 0.3451, (3.10a) 

[as{mr)]cv = 0.3182, (3.10b) 

4- Loop: 

0.2957 < as{mr) < 0.3477, (3.11a) 

[as{mr)]cv = 0.3200, (3.11b) 
Pade-Improved : 

0.2963 < as{mr) < 0.3489, (3.12a) 

[asimr)]cv = 0.3208. (3.12b) 

The ranges listed above progressively overlap the low end of the present experimental range 0.337 < 
as{mr) < 0.403. The central values (cv) displayed above are evolved down from qi^^Mz) = 0.118 with a 
five-flavour threshold at nth = 4.3 GeV. The lower bounds evolve from ^^(M^) = 0.115 with /ith = 4.0 
GeV, and the upper bounds evolve from as{Mz) = 0.121 with /i^h = 4.6 GeV. 

There are also some unexpected theoretical consequences arising from estimating the summation of 
the full /3-function series through use of an appropriately chosen Pade approximant [1] . The most general 
[2|2] approximant must yield 1 + Rix + R2x'^ + R^x^ + R4X^ as the first five terms of its Maclaurin 
expansion. For optimal generality, we require that Ri, R2, and R3 be given by (3.1b) and (3.2), but allow 
i?4 to be arbitrary. We then find that 



P{x) = -(3ox 



1 + aix + a2X 



(3.13) 



1 + bix + b2x'^ 
with {ai, 02) &i) ^2} linear in i?4 as follows: 
n/ = 3: 

ai = 7.1945 - 0.10261i?4, (3.14a) 

02 = -11.329 + 0.075643i?4, (3.14b) 

61 = 5.4168 - 0.10261i?4, (3.14c) 

62 = -25.430 + 0.25806i?4; (3.14d) 

n/ = 4: 

ai = 4.8401 - 0.11068i?4, (3.15a) 
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02 = -8.1842 + 0.10836ii:4, (3.15b) 

bi = 3.3001 - 0.11068i?4, (3.15c) 

62 = -16.314 + 0.21904i?4 (3.15d) 

Uf = 5: 

ai = 2.6793 - 0.12329i?4, (3.16a) 

02 = -6.19671 - 0.011245i24, (3.16b) 

bi = 1.4184 - 0.12329/^4, (3.16c) 

62 = -9.4599 + 0.14421i?4; (3.16d) 

Hf = 6: 

01 = 0.61085 - 0.18424i24, (3.17a) 

02 = -6.6275 - 0.24181i?4, (3.17b) 

61 = -0.31772 - 0.18424i?4, (3.17c) 

62 = -6.0423 - 0.057574i?4. (3.17d) 



For all Uf values listed above, the first positive zero of 1 + aix + a2x'^ in (3.13) is found to be above 
the first positive zero of l + feix + 62.T^. regardless of the choice for i?4. Consequently, the smallest positive 
zero of (3(x), if given by (3.13), is necessarily an ultraviolet fixed point and not an infrared fixed point, 
an inescapable result of the denominator sign change for x between and the first positive zero of (3.13). 
Moreover, for those values of R4 for which a second positive zero of 1 + aix + 02^^ is possible, we find 
from (3.14-17) that a second positive zero of 1 + bix + b2x'^ will also occur at some value of x between the 
two positive zeros of 1 + aix + a2X^. This ensures that neither positive zero of 1 + aix + 02^^ corresponds 
to an infrared fixed point. 

In Figure 2, a schematic diagram is presented showing different branches for the evolution of x(/x) 
anticipated from a /3-function (3.13) with the above-described alternation of positive denominator and 
numerator zeros, with the smallest positive zero occurring for the denominator. Zeros of 1 -f bix + b2X^ 
represent values of x for which x{^) has infinite slope. Zeros of 1 -|- aix + a2X^ are fixed-point values of 
X. As is evident from the figure, all such fixed points are necessarily ultraviolet, and the infrared region 
is inaccessible for values of fj, less than those corresponding to zeros of the denominator. Such behaviour 
suggests... 

1) ...the possible existence of a strong phase of QCD at short-distances, reflective of a nonzero ultraviolet 
fixed point, and 
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2) ... the inapplicability of a perturbative theory of quarks and gluons to the infrared region, specifically 
the region excluded from the domain of the first positive branch of x{fi) [Fig 2]. 

The former statement above may have ramifications for scenarios of dynamical electroweak symmetry 
breaking that usually involve a distinct technicolour group. The latter statement parallels old infrared 
slavery ideas, except that the inapplicability of perturbation theory to low energies is not seen to follow 
from the coupling constant growing infinite (or nonperturbatively large), as in infrared slavery, but from 
an explicit decoupling of the infrared region from the ultraviolet by virtue of /^-function singularities 
alternating with /3- function zeros (Fig. 2). We have verified that this alternation occurs in the most 
general [2|2] Pade-summation of the /3-function even when Uf = 0. 



4. Pade-Improvement of the Running Mass 

The running mass m(/i) satisfies the differential equation 

dm -f{x) 

where x = as{fJ,)/'7T, as in the previous section, and where f3{x) (= ji^dx/dii^) is given by (3.1a) and (3.2). 
The QCD MS anomalous mass dimension function 7(x) has been calculated to four-loop order [4,6]: 

^{x) = -x[l + Y,liA (4.2a) 

71 = 4.20833 - 0.138889n/ (4.2b) 

72 = 19.5156 - 2.28412n/ - 0.0270062n^ (4.2c) 

73 = 98.9434 - 19.1075^/ + 0.276163n^ - 0.00579322n^. (4.2d) 

Fade improvement of the square-bracketed expression within (4.2a) is straightforward via the methods of 
Section 2 by identifying R\,R2 and in (2.1) with 71,72 and 73. Using (2.12) we obtain the following 
APAP estimates for 74: 

(4.3) 
(4.4) 
(4.5) 
(4.6) 

One can obtain a solution for m[x(/i)] which includes the 0{x'^) Fade improvement of f3{x) and ^{x) by 
expressing 'y{x)/P{x) as a Maclaurin series in x, using APAP estimates(3.3-6) for P4 and (4.3-6) for 74. 
By truncating this series after x^, the differential equation (4.1) can be approximated by 

— ^ = \Po^ + dix + d2X^ + d3X^ + d^x'^] , (4.7) 
m dx \- J 



nf 


= 3 : 


74 = 


162.987 


Uf 


= 4 : 


74 = 


75.2349 


Uf 


= 5 : 


74 = 


12.5550 


Uf 


= 6 : 


74 = 


12.1820 
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with given by (3.2a), and with di given as follows: 





= 3 : 


di 


= 0.895063, d2 


= 1.94172, ds = 


2.88956, di = 


417.493; 


(4.8) 


Uf 


= 4 : 


di 


= 1.014131, d2 


= 1.74994, ds = 


0.0880435, di 


= -3.93256; 


(4.9) 


nf 


= 5 : 


di 


= 1.17549, d2 = 


= 0.809817, ds = 


-1.05016, di 


= -10.0138; 


(4.10) 


Uf 


= 6 : 


di 


= 1.39796, d2 = 


= 1.63266, d3 = - 


-6.84005, di = 


-- 142.769; 


(4.11) 



The values for di,d2 and ^3 are exactly determined by four-loop calculations of P{x) and 7(0;). The value 
of di is underlined to emphasize that it is determined from APAP estimates. The large values for di 
when Hf = 3 and Uf = 6 reflect correspondingly large values for ^i that are discussed in the previous 
section. 

The solution to (4.7) can be expressed in terms of x{fi) evaluated at two different values of /u: x{fii) = 
xi, x{n2) = X2, where x(/x) is the running coupling whose evaluation is discussed in the previous section. 
This solution to (4.7) is [4] 



m{x2) = m{xi)c{x2)/c{xi), 



(4.12a) 



where 



c{x) 



+ 
+ 



X 



l//3( 



{l + diX+[{dl + d2)/2 



{dl + 2d3 + 3did2)/6 



{df + 3dl + 6di + 6d?d2 + 8did3)/24] x^} 



(4.12b) 



Coefficients of x, x'^, and x^ are determined in full by known coefficients in the four- loop /3 and 7 
functions. The x'^ term is the lowest-degree term sensitive to Pade driven estimates of (3i and 74. We 
find the following set of expressions for c{x) from the di in (4.8-11); 



Uf = 3 : c{x) 



.4/9 



1 + 0.895063X + 1.37143a;^ 



+ 1.95168a:^ + 106.122a;^ 



(4.13) 



Uf = A : c{x) = 



,12/25 



1 + 1.01413X + 1.38920x^ 



+ 1.09052x''^ -00765827x1 



(4.14) 



nf = 5: c{x) = x^^/^^ [l -M.17549x -M.50071x2 
+ 0.172486x^ - 10.2813x^^1 , 



(4.15) 
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n/ = 6 : c{x) = x^l'' [l + 1.39796x + X.l^Mlx^ 
- 0.683486x3 + 33.7949x^ . 



(4.16) 



These same expressions are obtained to 0{x^) in ref. [4]; the effects of Pade- improvement reside entirely 
in the terms. 

It is important to recognise that these results ultimately derive from applying the asymptotic error 
formula (2.3) to the perturbative field-theoretical calculation of fi{x) and 7(0;), as argued in [2] and 
[3]. As in the previous section, the results (4.3-6) differ from those one would obtain using the fits of 
ref. [3] to the coefficients of within 74, particularly as 74 so extracted involves the near cancellation 
of large terms from successive values of k: from Table X of ref. [3] one finds for = 5 that 74 = 
530 - (143) • 5 + (6.67) • 5^ + (0.037) • 5^ - (8.54 • 10"^) • 5^ = -13.7. Small variations in these Pade- 
estimated coefficients can easily lead to positive values comparable to (4.5). 

It is also important to note that the application of the APAP algorithm at the field-theoretical level — 
i.e., to /?(x) and 7(0;) — is not equivalent to applying it to "perturbative" expressions which arc obtained 
by integrating over these functions. One could question, for example, whether the terms appearing in 
(4.13-16) might be obtainable by direct application of the APAP algorithm (2.12) to lower-degree terms 
in X. If we apply (2.12) directly to (4.13-16) using the explicit coefficients of .x, and in order to 
estimate the coefficients of x^, the x^-cocfficicnts we obtain are very different from those listed. Instead 
we obtain respectively forn/ = {3, 4, 5, 6} : 2.683x^, 0.7426x^, 0.01839x^, and 0.2850x^. This discrepancy 
is indicative of the inapplicability of the error formula (2.3) to the series in (4.12b), assuming (2.3) is 
applicable to the perturbative field-theoretical series (3.1a) and (4.2a). Such applicability is suggested by 
the predictions of five-loop terms for 0{N) j3- and 7-functions already noted in Section 2. 

Figure 3a-d display the relative impact on running quark masses of higher order corrections both 
augmented and unaugmented by Pade-improvement. Given an initial value mb(4.3 GeV) = 4.3 GeV [8], 
we evolve mb{n) up to = 175 GcV. Fig 3a indicates the evolution obtained via (4.1) from three-loop 
/3- and 7-functions — i.e. from truncation of the series (3.1a) and (4.2a) after i = 2. Fig. 3b displays 
the relative effects of higher-order corrections augmented and unaugmented by Pade-improvement. We 
first consider the unaugmented 4-loop case. The upper curve in Fig 3b is the ratio of nif, obtained from 
4- loop (3- and 7-functions to m^, obtained from 3- loop /?- and 7-functions (Fig 3a). As f^ increases from 4 
GeV to 175 GeV, the relative decrease in m},{^) from use of four-loop information is seen to be less than 



Pade improvement of four-loop results is displayed in the lower curve of Fig 3b. This curve is the 
ratio of a fully Pade-improved estimate of m(,(/u) to the three-loop calculation of m(,(/Lt) displayed in Fig 
3a. The full Pade-improvement is obtained via an APAP-algorithm determination of /?4 and 74, which is 
then utilized to construct Pade-summation [2|2] approximants as estimates of the aggregate effect of all 
higher-order terms in (3 and 7. The [2|2] approximant for /3(x) is given by (3.7c); the [2|2] approximant 
within 



is obtained via (2.15-19) using the n/ = 5 values of 71, 72, and 73 as well as the APAP-algorithm value 
for 74 given in (4.5). The full [2|2] approximants are then used to evaluate c(x) in (4.12a): c(x) = 
exp [/ (7(x)//3(x)) dx]. It is evident from Fig 3b that Pade-improvcmcnt docs not significantly alter 
mi,{fi) beyond a correction comparable to five-loop expectations; the relative change from such Pade 



0.1%. 



7(x) = — X 



1 + 1.19485X + 1.02765x^ 
1 - 2.31904X + 1.75664x2 



(4.17) 
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improvement is respectively within 0.1% (Fig 3b) and 0.01% of the unimproved three- and four- loop 
results. 

Fig. 3c displays the corresponding evolution of the charmed quark mass mcip) for 1.3 GcV < /v, < 
20 GeV, as obtained from (5- and 7- functions evaluated to two-, three-, and four- loop order as well 
as from Fade-improvement of the four-loop /?- and 7-functions, both below and above the five-flavour 
threshold. The initial value is taken to be mc(1.3 GeV) = 1.3 GeV for all four curves, and the five-flavour 
threshold is assumed to occur at 4.3 GeV [8]. The "Fade" curve is obtained by direct substitution of 
appropriate [2|2] Fade-summation (5- and 7-functions into (4.1). At the five-flavour threshold, we utilize 
the threshold-matching constraint [9] 

m^^/itfe) = m^-')(Mtfe) [1 + 0.2060 (a^\Mt/.)/7r)' 

+ (1.8229 + 0.0247n/)(ai''^^(;Ut/,)/7r)3]"^ (4.18) 

with Uf = 5 and fiffi = 4.3 GeV to generate the above-threshold initial condition for rric- Thus, the above- 
threshold portion of the Fade curve in Fig. 3c is obtained from this initial condition via substitution 
of (3.7c) and (4.17) into the differential equation (4.1). The below threshold portion of the Fade curve 
is obtained from the initial condition mc(1.3 GeV) = 1.3 GeV via substitution of the Fade-summation 
four-flavour /3-function (3.7b) and four-flavour 7-function 

[1 - 0.4541X - 1.3454a:^] 

^^^'^ ~ ~'^[l-4.1069x + 3.7090x2] (4.19) 

into the differential equation (4.1). The constraint (4.18) is also utilized to generate the four-loop curve 

in Fig. 3c, and [when taken to OTdev-(^a^J^^\fith)^ ] to generate the three-loop curve. As in Fig. 1, curves 
of successive order appear to converge (at least qualitatively) to the Fade estimate, which is almost 
indistinguishable from the four-loop curve. 

Fade- improvement effects are somewhat larger for light quarks (u, d, s). The evolution of light quarks 
from an initial value (normalized to unity) at = 1 GeV is displayed in Fig. 3d. This latter set of curves 
is obtained via utilization of (4.18) at both n/ = 4 and n/ = 5 flavour thresholds, which are respectively 
taken to be at 1.3 GeV and 4.3 GeV [8]. Below the 1.3 GeV four-flavour threshold, the Fade curve is 
generated via (4.1) using the ra/ = 3 Fade-summation 7-function 

, , [1 - 1.2373X - 1.8485x2] 

7(x) = -X7 ^ (4.20) 

^ [1 - 5.0289a; 4.7993x2] ^ ^ 

and the nj = 3 /3-function (3.7a). Between the four- and five-flavour thresholds we utilize (4.19) and 
(3.7b), and above the five-fiavour threshold, we utilize (4.17) and (3.7c) as before. At ^ = 5 GeV, Fig. 3d 
shows that there is a 1% difference between running masses obtained via unimproved and Fade-improved 
four-loop /3- and 7-functions. Once again, however, the distance between curves of successive order 
decreases as the order increases, giving the appearance of convergence towards the Fade-improved curve. 



5. Application to Higgs Decays 

Although the Higgs particle has yet to be directly observed, expressions for its decay into either two- 
gluons [10] or a bb pair [11] have been worked out with precision in perturbation theory. In much the 
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same way, knowledge of the Z° decay widths, whose precise values are important for bounds on standard- 
model parameters, preceded the discovery of the Z° itself. In this section, we apply Padc-improvcment 
to the decay processes H two gluons and H —>■ bb, and examine whether such improvement leads to 
detectable changes from the calculated rates obtained without such Pade improvement. 

5.1 Higgs — > Two Gluons 

The decay rate H ^ gg has been calculated to three-loop order in perturbation theory [10]: 

+ (156.808 - 5.70833 ln{ml / Mjj)) xjj + 0{x%)\ , (5.1) 

where xh = x{Mh) = C(s(Mh)/tt, and where Mh is assumed to be less than mt. Fade-improvement can 
enter this expression both in the actual value of as{MH) evolving from a Pade-improved /3-function, as 
well as in a Fade-driven estimate of the 0{x%) contribution to the square bracketed expression in (5.1). 
One cannot apply the APAP-algorithm of Section 2 to estimate this term because for a given value of 
Mh, only the coefficient of x (Ri = 17.9167) and [R2 = 156.808 - 5.70833 ln{m'^/M^)] arc known; 
the coefficient R3 of x^ is not known. One way to estimate R3 is to express 1 + Rix + R2X^ as a [1|1] Pade 
approximant, which upon expansion yields R^ = R^/ R\. A refinement on this estimate that is actually 
utilized to approximate /^s in ref. [2] is to assume in the asymptotic error formula (2.3) that a + 6 is small 
compared to A^-|- 1. One can then argue from (2.5) and (2.7) that 6^, — ^2/2, where 82 = (i?i — R2)/R2 is 
determined in full by Ri and i?2- Eq. (2.7) can then be rearranged to yield the following estimate of R^,: 

Rl/Rl 2Rl (r r,^ 

- TT^ - Rl + RiR2- ^^-^^ 

For different Mh values we plot in Fig. 4 the ratio of the Fade-improved H ^ gg rate to the rate 
obtained directly from (5.1), using the three-loop (3 function to obtain xh from the initial condition 
x{Mz) = 0.118/7r [8]. The Fade-improvement of the H ^ gg rate is obtained for a given choice of Mh 
first by evolving xh from the same initial condition via the [2|2] approximant (3.7c) for the /3-function, 
then by using (5.2) to estimate the 0{xff) term in (5.1), and finally by replacing the now-known cubic 
1 + RiXH + R2X% + R3XH ill (5.1) with its appropriate [2|1] Fade summation: 

I^R 3 1 + (i^l - R3/R2)X + {R2 - R3Rl/R2)x^ 

1 + Rix + R2X^ + Rsx-" ^ — . (5.3) 

1 - {R3/R2)X 

Figure 4 shows that such Fade improvement yields a 2.5% - 3% increase in the H ^ gg rate, with very 
little sensitivity to the Higgs mass. Such improvement is best understood to be an estimate of (unknown) 
higher order corrections to (5.1) that should be eventually testable against both experimental and future 
higher-order calculations of the H ^ gg rate. 

5.2 Higgs — > bb 

The decay rate Higgs — > bb has been calculated to four- loop order in perturbation theory [11]: 
r{H bb) = [sCpMHrnliMH)/ {^ttV2)^ 
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X { [l + (17/3)xif + 29.U67xjj + 41.7581x^] 

- i6ml{MH)/Mfj)[l + (20/3)xi^ + 14.62x|^]} (5.4) 
Full Pade-improvement of this expression for a given value of Mh (with Mh < mt) entails... 

1) ...determination of x{Mh) through use of (3.7), the [2|2] Pade-summation of the /^-function, to evolve 

x(/x) from an appropriate initial condition; e.g. ^^(Mj;) = 0.118 [8]; 

2) ...determination of mbiMn) through substitution into (4.1) of (3.7c) and (4.17), the [2|2] Pade 

summations for (3{x) and j{x), so as to evolve mh{fj,) from an appropriate initial condition; e.g. 
mb{A.3GeV) =4.3 GeV [8]; 

3) ...Pade-improvement and [2|2] Pade summation of the cubic expression in (5.4), 

' . im>. . ...14..^ . . ,5.) 

where (2.12) is used to generate an estimate of the coefficient [67.2472], and where (2.15 - 19) 
are used to generate the [2|2] approximant in (5.5), and 

4) ...Pade-improvement and [2|1] Pade summation of the quadratic expression in (5.4), 

where (5.2) is used to generate an estimate of the x^ coefficient [15.87], and where (5.3) is used to 
generate the [2|1] approximant in (5.6). 



The relative size of all these corrections, referenced to the rate calculated to the next-to-highest- 
known [three-loop] order in perturbation theory, is displayed in Fig. 5. The top curve is the ratio 
[r{H bb)U.ioop/[r{H ^ bb)]3 function of Mh- The 4-loop rate is obtained directly from 

(5.4), with mi,{MH) and xh{= x{Mh)) obtained from (3- and 7-functions that are truncated to zero after 
their 13^ and 73 contributions. The 3-loop rate is obtained by truncating off the highest order terms in 
(5.4) — specifically the 0{x^) term on the left-hand side of (5.5) and the Oix"^) term on the left-hand 
side of (5.6) - and by obtaining xh and mh{MH) from /?- and 7-functions that are truncated to zero after 
P2 and 72 contributions. The top curve shows a change of 0.02% to 0.09% in going fom three to four-loop 
order. 

The bottom curve compares the ratio ofr{H ^ bb), obtained by full Pade improvement of [T{H 
bb)], as described above, to [T{H — > 66)]3-;oop- It is evident from the figure that the two ratios are 
within 0.0001 of each other for all values of Mh below m^. In other words, Pade-improvcmcnt reduces 
the 4-loop result for r(H bb) by at most 0.01%. The effect that is seen seems to derive wholly from 
the Pade improvement of mfo(/Li), as is evident from comparison of Fig 5 to Fig 3b. Of course, such close 
agreement between (5.4) and its fully Pade-improved version suggests that the expression (5.4) is more 
than adequate for future comparison to experiment. Thus the purpose of the analysis presented here is 
really to demonstrate the robustness of (5.4) against Pade estimates of higher order corrections. 

The small size of corrections past even the 3-loop order is partly consequence of the small size of 
x{Mh) characterizing Higgs decay rates. Pade corrections are of much more interest when the magnitude 
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of X is larger, suggesting their usefulness in assessing the perturbative content of low-energy QCD - i.e. 
QCD sum rules. In the section that follows we will address how Pade-improvement can be utilized to 
estimate substantial higher-order corrections to sum rules relevant to scalar- and pseudoscalar-meson 
static properties. 



6. Perturbative Content of Scalar/Pseudoscalar QCD Sum Rules 

The resonance content of finite-energy and Laplace IZk QCD sum rules [12,13] is obtained from 
integrals over the imaginary part of current correlation functions n(s, //^) in the subcontinuum region 
(s < So): 



1 u 
^k{so) = - ImIl{s,so)s''ds, 

TT Jo 

7^fe(T, so) = - r ImU{s, l/T))s^e-'^ds. 



(6.1) 



(6.2) 



We consider here the purely perturbative content of the correlation function for scalar currents, which is 
presently known to 4-loop order [11]: 



—Imll{s, fi^) 

TT 

87r2 r ^ V TT 



oq -|- ailn ( — 

+ 62 ilni^ 



+ 



+ 



+ 



TT 



60 -|- biln 



s 

2\ / /..2' 



TT 



as{ix) 



Co -I- Clin I ~ I + '^2 M"- I — 



R4 + ...\. 



+ C3 /n 



(6.3) 



The problem we will address in this section is the computation of R4, which is necessary for the deter- 
mination of 0{ag) contributions to J-^ and TZ^. For three flavours Chetyrkin [11] has found that 



Oo = 17/3, ai = 2, 

60 = 31.8640, 61 = 31.6667, 62 = 17/4, 
Co = 89.1564, ci = 297.596, C2 = 229/2, cg = 9.20833. 
If we define w = s/fp, we can estimate i?4[tf] directly by substituting 
Ri [w] = uq — ai In w 
R2[w] = bo — bi In w + b2{ln w)'^ 



(6.4a) 
(6.4b) 
(6.4c) 

(6.5a) 
(6.5b) 
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-^sM = cq — ci In w + C2{ln luf' — c^{ln luf 



(6.5c) 



directly into (2.12). This enables one to determine explicitly 0{a!^) corrections to the sum-rules (6.1) and 
(6.2), even though R4{w\ determined in this way is manifestly not a fourth-order polynomial in ln{w). In 
particular, we easily find the C(af) contribution to J^k{so) to be 



AJPb(so) = 




2R4[w] dw 



(2059.4), 



(6.6a) 



A.Fi(so) = 




(6.6b) 



- ^ (^) («^^-^^)- (^•^^) 

The integrals in (6.6) have been evaluated numerically. This approach, however, ignores the known 
structural dependence of R4 on the variable w, 

R^lw] = do — diln w + d2{ln w)"^ — d^{ln w)^ + di{ln w)'^, (6.7) 

which may be important when one integrates over the ij;-variablc, as in (6.6). The ©(af) corrections to 
the first three finite energy sum rules are easily determined in terms of the constants di by substitution 
of (6.7) into the integrand of (6.1): 

j.T,r ^ 4 I «£(4^V (. ^di 2d2 ^ 2d3 ^ 8dA 



(6.10) 
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We can use the Pade algorithm (2.12) to estimate the coefHcients dj. To do so, we let Ri, R2, and R3 
be given by (6.5) for five representative values of w between zero and one: w = 

When w = 1, corresponding to s = sq in the finite-energy sum rule integrand (6.1), we see from (6.5) 
that Ri = ao, R2 = bo, R3 = cq. Using the APAP-algorithm (2.12), we find that R4I] = 251.442 = do- 
When w = e-V4 (s = 0.779so), we find from (6.5) that Ri = 37/6, R2 = 40.0463, R3 = 170.8554. Using 
(2.12) and (6.7), respectively, we then find that 



i?4[e-V4] = 699.398 = do + di/4 + ^2/16 + ^3/64 + ^4/256. 



Similarly we find the following results when w 



-1/2 (, 



O.6O6S0) 



w 



w 



(s = 0.135so): 



1/2] = 1389.82 = do + di/2 + ^2/4 + 0^3/8 + ^4/16, 
3652.36 = do + di + d2 + d3 + d^, 



i?4[e 
i?4[e-^] 

i?4[e-2] = 12804.9 = do + 2di + 4^2 + 8^3 + 16^4- 



(6.11) 

0.368so), and 

(6.12) 
(6.13) 
(6.14) 



We solve the four linear equations (6.11-14) for the four unknowns ^1,^2,^3,^4 using the value already 
obtained for do{= 251.422), and we obtain di = 1357.84, ^2 = 1634.53, ^3 = 404.630, d^ = 3.9097. Sub- 
stitution of these numbers into (6.8-10) yields results remarkably close to those obtained in (6.6). These 

results are listed in the underlined highest-order terms given below for the Pade-improved perturbative 

1/2 

content of the first three finite energy sum rules (x = q;s(sq )/tt): 



^o{so) 



167r2 



1 + 6.66667a; + 49.8223^2 



+ 302.110x^ + 2057.0x'^ 



(6.15a) 



^i{so) = 



87r2 



1 + 6.33333a; + 43.3640^2 



+ 215.846x^ + UMAx"^ 



(6.15b) 



^2(80) = 



34 



1 + 6.16667a; + 40.3119a;2 



327r2 

+ 178.731x^ + 833.52 x'^" . 



(6.15c) 



Pade corrections to the Laplace sum rules (6.2) arc not listed, as they are complicated by the occurrence 
of two scale variables (sq and r). However, such corrections are straightforward to obtain via integration 
of the 0{ag) term of (6.3), which we have already obtained via APAP estimates of do-4- 



R4 



251.44 + 1357.8 In 



+ 1634.5 In 



+ 404.63 In 



+ 3.9097 In 



(6.16) 
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It is worth noting that the underhned R4 terms in (6.15) are not very different from those one would 
obtain using either the R'^/R2 estimate suggested by a [2|1] approximant, as discussed following (2.10), or 
the APAP algorithm (2.12) applied directly to the known (non-underlined) terms of (6.15). The increase 
in the size of coefficients with increasing powers of x suggests the utility of a [2|2] Pade summation for 
these three expressions as an improvement over truncating off what may be substantial 0{x^) corrections 
to (6.15). By applying (2.15-19) to the three equations (6.15), we obtain the following [2|2] approximants 
to the full perturbative content of the first three finite-energy sum rules: 



34_ 
167r2 



1 + 3.8073a; + 6.8124x2 
1 - 2.8594a; - 23M7x^ 



(6.17a) 



^i(so) 



„3 

87r2 



1-h 2.3917.x -Ml. 308.x2 
1 - 3.9416X - 7.0930x2 



(6.17b) 



•^2(50) 



327r2 



1 + 2.2174X + 12.791x2 
1 - 3.9492X - 3.1670x2 



(6.17c) 



1 /2 

Both (6.15a) and (6.17a) arc indicative of a need for ,Sq to be substantially larger than 1 GeV for 

1/2 

finite-energy sum rules to be useful in the scalar and pseudoscalar channels. If Sq = 1 GeV, we see 
from Fig 1 that x{lGeV) = 0.153. For this value of x, each successive term in the square brackets of 
(6.15a) is approximately unity, indicative of nonconvergence. This is reflected by a near-vanishing of the 
denominator of (6.17a), implying a divergent result for the summation of the full perturbative series. If 

So = 3.24 GeV^, we see from Fig. 1 that x(sy^) = 0.10. The truncated series (6.15a) is then seen to yield 
a value that is only 87% of that obtained via the Pade summation (6.17a), indicative of the magnitude of 
the higher order terms missing from (6.15a). Note that a choice for sq near or somewhat above 3 GeV^ 
is suggested by Laplace sum-rule fits in both the pseudoscalar [14] and scalar [15] resonance channels. 

The finite energy sum rule J^q{sq) provides an example of how it is not enough just to have precise 
higher-order results. Even though (6.15a) includes 4-loop effects as well as an APAP-algorithm estimate 
of 5-loop effects, the five terms listed demonstrate only sluggish convergence for a realistic choice of sq. 
There is found to be enough of a difference between the truncated series (6.15a) and its Pade-summation 
(6.17a) to suggest the advisability of using the latter. 



7. Summary 

Using a Pade- motivated algorithm (2.12), we have estimated in Section 3 the five-loop contributions to 
the /3-function for rif = {3, 4, 5, 6}, and we have compared the evolution of oisit^) from ^ = Mz obtained 
from two-loop, three-loop, four-loop, and Pade-summation estimates of the full /3-function. Low energy 
values of obtained from the four-loop /3-functions with quoted flavour thresholds and appropriate 
threshold matching conditions are within 1% of those obtained from the Pade-summation /3-functions, a 
small effect compared to the much larger sensitivity of Q;s(l GeV) and asirur) to present uncertainties in 
as{Mz) and c- and b-quark flavour thresholds. 

We concluded Section 3 by extracting the most general set of [2|2] Pade-summation estimates of 
3, 4, 5, and 6 flavour QCD /3-functions whose Maclaurin expansions yield known four-loop results for 
their first four terms. For positive values of ag, these Pade-summation estimates of the /^-function were 
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shown to alternate denominator and numerator zeros, regardless of the size of the (presently unknown) 
five-loop term serving as a free parameter in the approximants. Such alternation necessarily implies that 
all positive numerator zeros represent ultraviolet fixed points, behaviour which, if applicable to the true 
/3-function, would decouple the (suitably defined) infrared region from perturbative QCD. 

In Section 4, we applied Pade-improvement methods to the running quark mass by estimating five- 
loop contributions to the 7-functions for 3, 4, 5, and 6 quark flavours. We then extracted an estimate for 
the O(x^) contribution to closed-form expressions for mq[as(/u)/7r] that had been obtained earlier [4,6] to 
0{x^). We compared the evolution of 3-loop, 4-loop, and Pade-summation estimates of rribdj,), once again 
finding very little relative difference (0.01%) between Pade-summation and four-loop determinations of 
mb{fJ,) over the range fj, < rrit, given identical five- flavour-threshold initial conditions. Corresponding 
agreement was still seen to occur at the 1% level for light quarks. 

In Section 5 we applied the results of the previous two sections to higher-loop calculations of the Higgs 
decay ratio H ^ gg and H — > 66, rates which are sensitive to running couplings and running masses, as 
well as higher- loop corrections that are polynomial in asimn)- The calculated three- loop H ^ gg rate 
is shown to be within 3% of our Pade-improvement estimate, given identical choices for M//, ^(M^), and 
mft-threshold initial conditions. Similarly, the calculated four-loop — > 66 rate is seen to differ from full 
Pade-improvement by at most 0.01%. 

All the results summarized up until this point are indicative of close agreement between known 
perturbation theory and Pade-approximant improvements intended to take into account higher-order 
effects. Consequently, the theoretical uncertainties associated with the truncation of any such calculations 
at the three-or-four-loop order are shown to be small. In Section 6, we considered quantities known 
to four-loop order for which this is not the case, the purely-perturbative content of QCD sum rules 
in scalar /pseudoscalar-resonancc channels. We constructed a Pade-algorithm estimate of the purely- 
perturbative O(a^) contribution to the imaginary part of scalar/pseudoscalar correlation functions, and 
we obtained [2|2] Pade-summation estimates of the all-orders content of the first three finite energy 
sum rules. We found the overall convergence of the primary sum rule to be doubtful for values of the 
QCD continuum threshold near sq = 1 GeV^. Even for sq above 3 GeV^, we found a greater-than-10% 
discrepancy between Pade-summation and four-loop-ordcr contributions to this sum rule, suggesting the 
existence of substantial theoretical uncertainties from higher-than-four-loop contributions. 
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I3i via (2.12) 




/34 via [2] 


(74/70) via (2.12) 


(74/70)*™' 


N = 


1 


1421.7 


1424.3 


1432.3 


135.1 


150.76 


N = 


2 


1941.7 


1922.3 


1943.8 


168.4 


191.89 


N = 


3 


2555.9 


2499.3 


2540.3 


203.1 


236.94 


N = 


4 


3267.9 


3158.8 


3225.6 


239.2 


285.94 



Table 1: Comparison of 0{N) /3-function and 7- function coefficients obtained via Fade estimates to those 
obtained via exact calculation. The "/34 via (2.12)" column displays estimates obtained via the APAP 
algorithm (2.12), based on Po-3 listed in [7]. The column displays the results for obtained 

in [7] by explicit calculation. The "/34 via [2]" column lists the Fade estimates given in Table 3 of ref. 
[2], involving knowledge of the A^^ dependence and a fit of the overall A^-dependence of (34, as well as a 
simplified asymptotic error formula. The "(74/70) via (2.12)" column displays estimates of 74/70 in the 
anomalous mass dimension 7-function obtained via the APAP algorithm (2.12), based on prior coefficients 
listed in [7]. The final column, "(74/70)*™'" are coefficients explicitly calculated in [7]. 

Figure Captions: 

Figure 1: A comparison of the evolution of x{p) = as{p,)/TT for rrit > fJ, > 1 GeV obtained from the 
evolution of two-loop (2L), three-loop (3L), four-loop (4L), and [2|2] Pade-summation /3-functions. 
All curves are generated from the initial condition as(Mz) = 0.118, with five- and four-fiavour 
threshold matchings occurring at 4.3 GeV and 1.3 GeV, respectively. 

Figure 2: Schematic behaviour of x(/x) obtained from a [2|2] Pade-summation estimate of the /3-function 
whose positive numerator and denominator zeros alternate. The denominator zeros are denoted by 

Xdi and Xd2, the numerator zeros arc denoted by x„i and a;„2; and the alternation of zeros is 
consistent with the smallest positive zero being a zero of the denominator: < Xdi < Xni < X(i2 < 
Xn2- The value /idi is defined such that x{iidi) = Xdi; values oi jj, < fj,di are outside the domain of 

Figure 3a: Evolution of the running mass mb(/n) obtained from three-loop (3L) /3- and 7-functions from 
the initial condition m6(4.3 GeV) = 4.3 GeV. 

Figure 3b: A comparison of the ratio of mb{ii) obtained from four-loop (4L) and [2|2] Pade-summation 
13- and 7-functions to mh{n) obtained from three-loop (3L) /?- and 7-functions, as shown in Fig. 3a, 
given the initial condition m6(4.3 GeV) = 4.3 GeV. 

Figure 3c: A comparison of the evolution of mc{^i) obtained from evolution of two-loop (2L), three-loop 
(3L) four-loop (4L) and [2|2] Pade-summation (3- and 7-functions. 

Figure 3d: Masses of light (u, d, s) quarks from 2-loop, 3-loop, 4-loop, and [2|2] Pade-summation /?- 
and 7-functions. Masses at // = 1 GeV are normalized to unity. 

Figure 4: A comparison of the ratio of the Higgs decay rate into two gluons obtained via Fade- improvement 
discussed in the text to the three-loop (3L) expression for the same rate. 

Figure 5: A comparison of the ratio of the Higgs decay rate into a bb pair obtained to four-loop (4L) order 
and through subsequent Pade-improvement (as described in the text) to the same rate obtained to 
three-loop (3L) order. 
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